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We consider the cosmological implications of light travelling faster in the early Universe. We propose 
a prescription for deriving corrections to the cosmological evolution equations while the speed of 
light c is changing. We then show how the horizon, flatness, and cosmological constant problems 
may be solved. We also study cosmological perturbations in this scenario and show how one may 
solve the homogeneity and isotropy problems. As it stands, our scenario appears to most easily 
produce extreme homogeneity, requiring structure to be produced in the Standard Big Bang epoch. 
Producing significant perturbations during the earlier epoch would require a rather careful design 
of the function c(t). The large entropy inside the horizon nowadays can also be accounted for in 
this scenario. 
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PACS Numbers: 98.80. Cq, 98.80.-k, 95.30.Sf 
I. THE PUZZLES OF THE BIG BANG MODEL 

Cosmologists have long been dissatisfied with the 
"Standard Big Bang" (SBB) model of the Universe. This 
is not due to any conflict between the big bang theory 
and observations, but because of the limited scope of- 
fered by the SBB to explain certain striking features of 
the Universe. From the SBB perspective the homogene- 
ity, isotropy, and "flatness" of the Universe, and the pri- 
mordial seeds of galaxies and other structure are all fea- 
tures which are "built in" from the beginning as initial 
conditions. Cosmologists would like to explain these fea- 
tures as being the result of calculable physical processes. 
A great attraction of the Inflationary Cosmologies Q is 
that they address these issues by showing on the basis of 
concrete calculations that a wide variety of initial condi- 
tions evolve, during a period of cosmic inflation, to re- 
flect the homogeneity, isotropy, flatness and perturbation 
spectrum that we observe today. 

So far, all attempts to achieve this kind of improve- 
ment over the SBB have wound up taking the basic infla- 
tionary form, where the observable Universe experiences 
a period of "superluminal" expansion. This is accom- 
plished by modifying the matter content of the Universe 
in such a way that ordinary Einstein gravity becomes 
repulsive and drives inflationary expansion. This pro- 
cess is in many ways remarkably straightforward and has 
found numerous realizations over the years ( [^^), etc), 
although it might still be argued that a truly compelling 
microscopic foundation for inflation has yet to emerge. 

One interesting question is whether inflation is the 
right solution to the cosmological puzzles. Is inflation 
really what nature has chosen to do? When this mat- 
ter is discussed there is a notable absence of any real 
competition to inflation, and this must be counted in in- 
flation's favour. However, we believe the picture would 
become much clearer if some kind of debate along these 
lines were possible. To this end, we discuss here a possi- 
ble alternative to inflationary cosmology which, while not 
as well developed as today's inflationary models, might 
lead to some illuminating discussion. 

In this alternative picture, rather than changing the 



matter content of the Universe, we change the speed of 
light in the early Universe. We assume that the Uni- 
verse matter content is the same as in the SBB, that is, 
the Universe is radiation dominated at early times. We 
also assume that Einstein's gravity is left unchanged, in 
a sense made precise in Section [V. The geometry and 



expansion factor of the Universe are therefore the same 
as in the SBB. However the local speed of light, as mea- 
sured by free falling observers associated with the cosmic 
expansion, varies in time, decelerating from a very large 
value to its current value. 

We discuss below how Varying Speed of Light (VSL) 
models might resolve the same cosmological puzzles as 
inflation, and offer a resolution to the cosmological con- 
stant problem as well. We shall not dwell on the possible 
mechanisms by means of which the speed of light could 
have changed. Rather we wish to concentrate on the 
conditions one should impose on VSL models for their 
cosmological implications to be interesting. This phe- 
nomenological approach should be regarded as a curios- 
ity, which, we hope, will prompt further work towards an 
actual theory in which the physical basis of VSL models 
is realized. 

One may doubt that such a self-consistent theory could 
ever be constructed. We therefore feel forced to tran- 
scend the scope of this paper, and discuss essential as- 
pects of such a theory. We find it befitting to start our 
discussion with an assessment of the experimental mean- 
ing of a varying c (Section [nj). We also need to be more 
specific about VSL theories in order to tackle the flatness, 
cosmological constant, homogeneity, and entropy prob- 
lems. In Section [V we state what is actually required 
from any VSL theory to solve these problems. However in 
Appendix I we lay out the foundations for such a theory. 



II. THE MEANING OF A VARIABLE SPEED OF 
LIGHT 

We first address the question of the meaning of a vary- 
ing speed of light. Could such a phenomenon be proved 
or disproved by experiment? Physically it does not make 
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sense to talk about constancy or variability of any di- 
mensional "constant" . A measurement of a dimensional 
quantity must always represent its ratio to some stan- 
dard unit. For example, the length of my arm in meters 
is really the dimcnsionlcss quantity given by the ratio of 
the arm length to the length of a meter stick. If the ratio 
varied, one could interpret this as a variation in cither 
(or both) of the two lengths. In familiar situations, there 
is usually a preferred interpretation which distinguishes 
itself by giving a simpler view of the world. Choosing a 
given person's arm as a standard of length would require 
a whole range of simple objects to undergo peculiar dy- 
namics, whereas assuming the meter stick to be constant 
would usually give a much simpler picture. 

None the less, a given theory of the world requires di- 
mensional parameters. If these parameters varied, how 
would that process show up in experiments? Suppose we 
set out to measure the speed of light. For this one needs 
a length measure (rod) and a clock. In a world described 
by a theory with time varying dimensional parameters, 
it is quite possible that the rods and clocks, as well as 
the photon speeds, could all vary. Because measurements 
are fundamentally dimensionless, the experimental result 
will only measure some dimensionless combination of the 
fundamental constants. Let us sketch a simple illustra- 
tion: Suppose we measure time with an atomic clock. 
Taking the Rydberg energy (E R = m e e A /2(4ne ) 2 h 2 ) to 
represent the dependence of all atomic energy levels on 
the fundamental constants, the oscillation period of the 
atomic clock will be cx h/E R . Likewise, taking the Bohr 
radius (oq — Aire^h 2 /m e e 2 ) to reflect the relationship be- 
tween the lengths of ordinary objects (made of atoms) 
and the fundamental constants, the length of our rod is 
cx ao- Thus a measurement of c with our equipment is 
really a measurement of the dimensionless quantity 
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essentially the fine structure constant. We could of course 
use other equipment which depends in different ways on 
the fundamental dimensionless constants. For example, 
pendulum clocks will necessarily involve Newton's con- 
stant G. Different experiments will result, which measure 
different dimensionless combinations of the fundamental 
dimensional constants. Our conclusion that physical ex- 
periments are only sensitive to dimensionless combina- 
tions of dimensional constants is hardly a new one. This 
idea has been often stressed by Dicke (eg. and we 

believe this is not controversial. 

Thus, speaking in theoretical terms of time varying di- 
mensional constants can lead to problems. To give an 
historical example, papers |]l8| , |l9| were written claiming 
stringent experimental upper bounds on the time vari- 
ability of the dimensional quantity he. In these the prod- 
uct EX was found to be the same for light emitted at very 
different redshifts. From the deBroglie relation he = EX 
one infers the constancy of he. Bekenstein gives an illumi- 
nating discussion of the fallacy built into this argument 



p0[ . Built into E cx 1/a and A cx a is the assump- 
tion that he is constant, for otherwise the wavevector 
and the momentum vector p M could not both be parallel 
transported. Hence the experimental statement that he 
is constant is circular. 

What would we do therefore if we were to observe 
changing dimensionless quantities? Any theory explain- 
ing the phenomenon would necessarily have to make use 
of dimensional quantities. It would a priori be a matter of 
choice, prejudice, or convenience to decide which dimen- 
sional quantities are variable and which are constant (as 
we mentioned in the illustration above). There would be 
a kind of equivalence, or duality between theories based 
on any two choices as far as dimensionless observations 
are concerned. However, the equations for two theories 
which are observationally equivalent, but which have dif- 
ferent dimensional parameters varying, will in general not 
look the same, and again simplicity will end up being an 
important factor in making a choice between theories. In 
what follows, we will prefer to work with models which 
have the simplicity of "minimal coupling" . 

Let us illustrate this point with a topical example. 
There has been a recent claim [ p2| of experimental ev- 
idence for a time changing fine structure constant a = 
e 2 /(4:Trhc). Although the ongoing chase for systematics 
precludes any definitive conclusions, let us assume for the 
purpose of the argument that the effect is real. 

In building a theory which explains a variable a we 
must make a decision. We could postulate that electric 
charge changes in time, or, say, that he must change in 
time. Bekenstein g^j constructs a theory based on the 
first alternative. He postulates a Lorentz invariant ac- 
tion, which does not conserve electric charge. Our the- 
ory is based on the second choice. We postulate breaking 
Lorentz invariance, a changing he, and consequently non- 
conservation of energy. Any arguments against the ex- 
perimental meaning of a changing c can also be directed 
at Bekensteins' changing e theory, and such arguments 
are in both cases meaningless. In both cases the choice of 
a changing dimensional "constant" reverts to the postu- 
lates of the theory and is not, a priori, an experimental is- 
sue. The observables are always dimensionless. However, 
the minimally coupled theories based on either choice are 
not dual (as we shall point out in Appendix I) . For this 
reason one might prefer one formulation over the other. 

Finally, and on a different tone, suppose that future 
experiments were to confirm that not only a changes in 
time, but also that there are time variations in dimen- 
sionless coupling constants based on other interactions, 
en = g 2 /(hcj^\. Suppose further that the ratios between 
the various constants, ry = cti/ctj, were observed to be 



In writing these constants we have assumed that the cou- 
plings of these interactions are defined in terms of "charges" 
(with dimensions of [E] 1/2 [L] 1/2 ). 
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FIG. 1. Conformal diagram (light at 45°) showing the hori- 
zon structure in the SBB model. Our past light cone contains 
regions outside each others' horizon. 

constant. Choosing what dimensional constants were in- 
deed constants would still be a matter of taste. One could 
still define a theory in which the various charges gi change 
in time, with fixed ratios, and he remains constant. How- 
ever it would perhaps start to make more sense, merely 
for reasons of simplicity, to postulate instead a changing 

he. 

Therefore, even though a variable c cannot be made a 
dimensionless statement, evidence in favour of theoretical 
models with varying c could be accrued if the other 
changed, with fixed ratios. 

III. COSMOLOGICAL HORIZONS 

Perhaps the most puzzling feature of the SBB is the 
presence of cosmological horizons. At any given time 
any observer can only see a finite region of the Universe, 
with comoving radius rt = cq, where rj denotes confor- 
mal time, and c the speed of light. Since the horizon size 
increases with time we can now observe many regions in 
our past light cone which are causally disconnected, that 
is, outside each others' horizon (see Fig. |l|). The fact 
that these regions have the same properties (eg. Cosmic 
Microwave background temperatures equal to a few parts 
in 10 5 ) is puzzling as they have not been in physical con- 
tact. This is a mystery one may simply relegate to the 
setting up of initial conditions in our Universe. 

One may however try to explain these very peculiar 
initial conditions. The horizon problem is solved by in- 
flationary scenarios by postulating a period of accelerated 
or superluminal expansion, that is, if a is the expansion 
factor of the Universe, a period with a > 0. The Fried- 
man equations require that the strong energy condition 
p + 3p/c 2 > must then be violated, where pc 2 and p 
are the energy density and pressure of the cosmic matter. 
This violation is achieved by the inflaton field. If a > 
for a sufficiently long period one can show that cosmo- 
logical horizons are a post-inflation illusion, and that the 
whole observed Universe has in fact been in causal con- 
tact since an early time. 

A more minimalistic way of solving this problem is to 
postulate that light travelled faster in the Early Universe. 
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FIG. 2. Diagram showing the horizon structure in a SBB 
model in which at time t c the speed of light changed from c~ 
to c + <C c~ . Light travels at 45° after t c but it travels at 
a much smaller angle with the space axis before t c . Hence 
it is possible for the horizon at t c to be much larger than 
the portion of the Universe at t c intersecting our past light 
cone. All regions in our past have then always been in causal 
contact. 

Suppose there was a "phase transition" at time t c when 
the speed of light changed from c~ to c + . Our past light 
cone intersects t — t c at a sphere with comoving radius 
r = c + (?;o — r) c ), where r]o and r\ c are the conformal times 
now and at t c . This is as much of the Universe after the 
phase transition as we can see today Q. On the other 
hand the horizon size at t c has comoving radius = 
c~r\ c . If c~/c + ^> t]o/t] c , then r -C r/,, meaning that the 
whole observable Universe today has in fact always been 
in causal contact (see Fig. ^) . Some simple manipulations 
show that this requires 

c~ 1 1 T + 

lo Sio ^jr > 32 - - log 10 z eq + - log 10 (2) 

where z eq is the redshift at matter radiation equality, and 
T c + and Tp are the Universe and the Planck tempera- 
tures after the phase transition. If T c + w Tp this implies 
light travelling more than 30 orders of magnitude faster 
before the phase transition. It is tempting, for symmetry 
reasons, simply to postulate that c~ = oo but this is not 
strictly necessary. 

IV. A PRESCRIPTION FOR MODIFYING 
PHYSICAL LAWS WHILE THE SPEED OF 
LIGHT IS VARYING 

Hidden in the above argument is the assumption that 
the geometry of the Universe is not affected by a chang- 
ing c. We have allowed a changing c to do the job nor- 
mally done by "superluminal expansion" . To enhance 
this effect we have forced the geometry to still be the 
SBB geometry We now elaborate on this assumption. 
We will propose a prescription for how, in general, to 
modify gravitational laws while c is changing. This pre- 
scription is merely the one we found the most fertile. In 
Appendix I we describe in detail a theory which realizes 
this prescription. 
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The basic assumption is that a variable c does not in- 
duce corrections to curvature in the cosmological frame, 
and that Einstein's equations, relating curvature to stress 
energy, are still valid. The rationale behind this postulate 
is that c changes in the local Lorentzian frames associ- 
ated with cosmological expansion. The effect is a special 
relativistic effect, not a gravitational effect. Therefore 
curvature should not feel a changing c. 

The previous statement is not covariant. However in- 
troducing a function c(t) is not even Lorentz invariant. 
So it is not surprising that a favoured gauge, or coor- 
dinate choice, must be made, where the function c(t) is 
specified, and in which the above postulate holds true. 
The cosmological frame (with the cosmological time t) 
provides such a preferred frame. 

In a cosmological setting the postulate proposed im- 
plies that Friedman equations remain valid even when 
c ^ 0: 
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(3) 
(4) 



where, we recall, pc 2 and p are the energy and pressure 
densities, K — 0, ±1 and G the curvature and the gravi- 
tational constants, and the dot denotes a derivative with 
respect to proper time. If the Universe is radiation dom- 
inated, p — pc 2 /3, and we have as usual a oc t 1 / 2 . We 
have assumed that a frame exists where c = c(t), and 
identified this frame with the cosmological frame. 

The assumption that Einstein's equations remain unaf- 
fected by decelerating light carries with it an important 
consequence. Bianchi identities apply to curvature, as 
a geometrical identity. These then imply stress energy 
conservation as an integrability condition for Einstein's 
equations. If c ^ 0, however, this integrability condition 
is not stress energy conservation. Source terms, propor- 
tional to c/c, come about in the conservation equations. 

Seen in another way, the conservation equations im- 
ply an equation of motion for free falling point particles. 
This is normally the geodesic equation, but now source 
terms will appear in the geodesic equation. Clearly a vi- 
olation of the weak equivalence principle is implied while 
c is changing §. This, of course, does not conflict with 
experiment, as we take c ^ only in the Early Universe, 
possibly for only a very short time (such as a phase tran- 
sition). 

Although this is a general remark we shall be concerned 
mostly with violations of energy conservation in a cosmo- 
logical setting. Friedman equations can be combined into 
a "conservation equation" with source terms in c/c and 
G/G: 
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In a flat Universe (K — 0) a changing c does not vio- 
late mass conservation. Energy, on the other hand, is 



proportional to c 2 . If, however, K ^ not even mass is 
conserved. 

In Eqn. | we have included the effects of G under the 
same postulate merely for completness. In such a formu- 
lation VSL does not reduce to Brans Dicke theory when 
c = 0, and G ^ 0. This is because we postulate that 
Fricdmann equations remain unchanged, which implies 
that the conservation equations acquire terms in c and 
G. In Brans Dicke theory one postulates exactly the op- 
posite: the conservation equations must still be valid, so 
that the weak equivalence principle is satisfied. While we 
could have taken this stance for c as well we feel that vio- 
lation of energy conservation is the hallmark of changing 
c. Variable c must break Poincare invariance, for which 
energy is the Noether current. Barrow has proposed 
a formulation of VSL which has the correct Brans Dicke 
limit. 



V. THE FLATNESS PUZZLE 

We now turn to the flatness puzzle. The flatness puzzle 
can be illustrated as follows. Let p c be the critical density 
of the Universe: 



Pc = 



8irG \ a 



(6) 



that is, the mass density corresponding to K = for a 
given value oid/a. Let us define e = 17 — I with £1 = p/p c . 
Then 



\P Pc 



(7) 



If p = wpc 2 (with w = 0), using Eqns.(||), (§), and (||) 
we have: 



and so 



P „a / n G „c e 
£ = -3~ (1+w -- + 2~— - 
p a G cl + e 

^ = --(2+(l + e)(l + 3«,))-£ 
p c a G 



e = (l + e)e-(l + 3iu) + 2-e 
a c 



(8) 
(9) 

(10) 



In the SBB e grows like a 2 in the radiation era, and like a 
in the matter era, leading to a total growth by 32 orders 
of magnitude since the Planck epoch. The observational 
fact that e can at most be of order I nowadays requires 
that either e = strictly, or an amazing fine tuning must 
have existed in the initial conditions (e < I0~ 32 at t — 
tp). This is the flatness puzzle. 

The e = solution is in fact unstable for any matter 
field satisfying the strong energy condition I + 3u> > 0. 
Inflation solves the flatness problem with an inflaton field 
which satisfies I + 3u> < 0. For such a field e is driven 



4 



towards zero instead of away from it. Thus inflation can 
solve the flatness puzzle. 

As Eqn. [To| shows a decreasing speed of light (c/c < 0) 
would also drive e to 0. If the speed of light changes in a 
sharp phase transition, with \c/c\ 3> a/a, we can neglect 



S = 



the expansion terms in Eqn. 10, Then e/e = 2c/ c so that 



e oc c . A short calculation shows that the condition (|2j) 
also ensures that e < 1 nowadays, if e w 1 before the 
transition. 

The instability of the K ^ Universes while c/c < 
can be expected simply from inspection of the non con- 
servation equation Eq. (|^). Indeed if p is above its criti- 
cal value, then K = 1, and Eq. (||) tells us that mass is 
taken out of the Universe. If p < p c , then K = — 1, and 
then mass is produced. Either way the mass density is 
pushed towards its critical value p c . In contrast with the 
Big Bang model, during a period with c/c < only the 
K = Universe is stable. 

Note that with the set of assumptions we have used a 
changing G cannot solve the flatness problem (cf . p|-[l0t ) . 

We have assumed in the previous discussion that we 
are close, but not fine-tuned, to flatness before the tran- 
sition. It is curious to note that this need not be the 
case. Suppose instead that the Universe acquires "natu- 
ral initial conditions" (eg. e w 1) well before the phase 
transition occurs. If such Universes are closed they rec- 
ollapse before the transition. If they are open, then they 
approach e = — 1. This is the Milne Universe, which is 
our case (constant G) may be seen as Minkowski space- 
time. Such a curvature dominated Universe is essentially 
empty, and a coordinate transformation can transform it 
into Minkowski space-time. Inflation cannot save these 
empty Universes, as can be seen from Eqn. [H^ Indeed 
even if 1 + 3w < the first term will be negligible if 
e w — 1. This is not true for VSL: the second term will 
still push an e = — 1 Universe towards e = 0. 

Heuristically this results from the fact that the viola- 
tions of energy conservation responsible for pushing the 
Universe towards flatness do not depend on there being 
any matter in the Universe. This can be seen from in- 
spection of Eqn. (||). 

In this type of scenario it does not matter how far be- 
fore the transition the "initial conditions" are imposed. 
We end up with a chaotic scenario in which Darwinian 
selection gets rid of all the closed Universes. The open 
Universes become empty and cold. In the winter of these 
Universes a phase transition in c occurs, producing mat- 
ter, and leaving the Universe very fine tuned, indeed as 
an Einstein deSitter Universe (EDSU). 



VI. THE COSMOLOGICAL CONSTANT 
PROBLEM 

There are two types of cosmological constant problems, 
and we wish to start our discussion by differentiating 
them. Let us write the action as: 



1 {R + 2A 1 ) 
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where Cm is the matter fields Lagragian. The term in Ai 
is a geometrical cosmological constant, as first introduced 
by Einstein. The term in A 2 represents the vacuum en- 
ergy density of the quantum fields [jll). Both tend to 
dominate the energy density of the Universe, leading to 
the so-called cosmological constant problem. However 
they represent two rather different problems. We shall 
attempt to solve the problem associated with the first, 
not the second, term. Ususally one hopes that the sec- 
ond term will be cancelled by an additional couter-term 
in the Lagrangian. In the rest of this paper it is the 
geometrical cosmological constant that is under scrutiny. 

If the cosmological constant A / then the argument 
in the previous section still applies, with p = p m + p\, 
where p m is the mass density in normal matter, and 
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is the mass density in the cosmological constant. One 
still predicts £l m + ^A = 1, with Q m — p m /p c and Q\ 
f> \ /),.. However now we also have 
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If A is indeed a constant then from Eq. (|12| 
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If we define e\ = pa/ Pm we then find, after some straight- 
forward algebra, that 



, a> , * c 1 + 6a 
e A = £A 3- 1 + w + 2— — - 
' a c 1 + e 



(15) 



Thus, in the SBB model, e\ increases like a 4 in the ra- 
diation era, like o 3 in the matter era, leading to a total 
growth by 64 orders of magnitude since the Planck epoch. 
Again it is puzzling that ea is observationally known to 
be at most of order 1 nowadays. We have to face another 
fine tuning problem in the SBB model: the cosmological 
constant problem. 

If c = the solution ea = is in fact unstable for any 
w > — 1. Hence violating the strong energy condition 
1 + 3w > would not solve this problem. Even in the 
limiting case w — —1 the solution e\ = is not an attrac- 
tor: 6a would merely remain constant during inflation, 
then starting to grow like a 4 after inflation. Therefore 
inflation cannot "explain" the small value of ea , as it can 
with e, unless one violates the dominant energy condition 
w > -1. 

However, as Eqn. ( p"5| ) shows, a period with c/c <C 
would drive e\ to zero. If the speed of light changes 
suddenly (|c/c| 3> a/a) then we can neglect terms in 
a/ a, and so 
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e A (l + £A) cl + e 
which when combined with e/e = 2c/c leads to 
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(16) 



(17) 



The exact constraint on the required change in c depends 
on the initial conditions in e and eA- In any case once 
both e w 1 and ea ~ 1 we have £a oc c 2 . Then we 
can solve the cosmological constant problem in a sudden 
phase transition in which 
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logic — > 64 - - log 
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(18) 



This condition is considerably more restrictive than (0), 
and means a change in c by more than 60 orders of mag- 
nitude, if T c + Ri Tp . Note that once again a period with 
G/G would not solve the cosmological constant problem. 

Equations (|l0| ) and ( |l5|) are the equations one should 
integrate to find conditions for solving the flatness and 
cosmological constant problems for arbitrary initial con- 
ditions and with arbitrary curves c(t). They generalize 
the conditions §) and (|l|) which are valid only for a 
starting point with e w 1 and ea ~ 1 and for a step 
function c{t). 

As in the case of the flatness problem we do not need 
to impose "natural initial conditions" (ea ~ 1) just be- 
fore the transition. These could have existed any time 
before the transition, and the argument would still go 
through, albeit with a rather different overall picture for 
the history of the Universe. 

If £A ~ 1 well before the transition, then the Uni- 
verse soon becomes dominated by the cosmological con- 
stant. We have inflation! The curvature and matter will 
be inflated away. We end up in a de-Sitter Universe. 
When the transition is about to occur it finds a flat Uni- 
verse (e = 0), with no matter (p m — 0), and with a 
cosmological constant. If we rewrite Eqn.(|l5|) in terms 
of e m = p m /pA, f° r e — and \c/c\ ^ a/a, we have 
e m = — 2(c/c)(l + e m ). Integrating leads to l + e m oc c~ 2 . 
We conclude that we do not need the presence of any 
matter in the Universe for a VSL transition to convert a 
cosmological constant dominated Universe into a EDSU 
Universe full of ordinary matter. This can be seen from 
Eqns. (^3|)-(|l4]). A sharp decline in c will always dis- 
charge any vacuum energy density into ordinary matter. 

We stress the curious point that in this type of scenario 
the flatness problem is not solved by VSL, but rather by 
the period of inflation preceding VSL. 



VII. THE HOMOGENEITY OF THE UNIVERSE 

Solving the horizon problem by no means guarantees 
solving the homogeneity problem, that is, the uncanny 
homogeneity of the currently observed Universe across 



many regions which have apparently been causally dis- 
connected. Although solving the horizon problem is a 
necessary condition for solving the homogeneity problem, 
in a generic inflationary model solving the first causes se- 
rious problems in solving the latter. Early causal contact 
between the entire observed Universe allows equilibration 
processes to homogenize the whole observed Universe. It 
is crucial to the inflation picture that before inflation the 
observable universe in well inside the Jeans length, and 
thus equilibrates toward a homogeneous state. However 
no such process is perfect, and small density fluctuations 
tend to be left outside the Hubble radius, once the Uni- 
verse resumes its standard Big Bang course. These fluc- 
tuations then grow like a 2 during the radiation era, like 
a during the matter era, usually entailing a very inhomo- 
geneous Universe nowadays. This is a common flaw in 
early inflationary models ||l2|| which requires additional 
fine-tuning to resolve. 

In order to approach this problem we study in Ap- 
pendix II the effects of a changing c on the theory of 
scalar cosmological perturbations |13|. The basic re- 
sult is that the comoving density contrast A and gauge- 
invariant velocity v are subject to the equations: 



A' 



a' 

3w— 4 


1) 


A = 


a 










)- 








kc 2 


w r 


— kc 


+ 1 + 


w 



(1 + w)kv - 2-wII T 
a 

' c 2 s k 3 a' fa' d_ 
k 1 + w 2k a \ a c 

' 2/3 3 



(19) 



1 + w k 2 c 2 



wU T (20) 



where k is the wave vector of the fluctuations, and T is 
the entropy production rate, Hp the anisotropic stress, 
and c s the speed of sound, according to definitions spelled 
out in Appendix II. 

In the case of a sudden phase transition Eqn. ( |l9| ) 
shows us that A oc c, regardless of the chosen equations 
of state for T and Ht- Hence 



AJ 
A" 



(21) 



meaning a suppression of any fluctuations before the 
phase transition by more than a factor of 10 -60 if con- 
dition (jl^) is satisfied. The suppression of fluctuations 
induced by a sudden phase transition in c can be intu- 
itively understood in the same fashion as the solution to 
the flatness problem. Mass conservation violation ensures 
that only a Universe at critical mass density is stable, if 
c/c <C 0. But this process occurs locally, so after the 
phase transition the Universe should be left at critical 
density locally. Hence the suppression of density fluctu- 
ations. 

We next need to know what are the initial conditions 
for A and v. Suppose that at some very early time ti one 
has c/c = and the whole observable Universe nowadays 
is inside the Jeans length: rj <C Cj?/j/v3. The latter 
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condition is enforced as a byproduct of solving the hori- 
zon problem. The whole observable Universe nowadays 
is then initially in a thermal state. What is more each 
portion of the Universe can be described by the canonical 
ensemble and so the Universe is homogeneous apart from 
thermal fluctuations . These are characterized by the 
mass fluctuation 



2 



(6M 2 



4fc b T, 
Mc? 



(22) 



Converted into a power spectrum for A this is a white 
noise spectrum with amplitude 



P A (fc)H|A(fc) 2 |)(x 



(23) 



What happens to a thermal distribution, its tempera- 
ture, and its fluctuations, while c is changing? In ther- 
mal equilibrium the distribution function of particle en- 
ergies is the Planck distribution P(E) = l/(e E / kbT - 1), 
where T is the temperature. When one integrates over 
the whole phase space, one obtains the bulk energy den- 
sity pc 2 oc (fcf,T) 4 / (he) 3 . Let us now consider the time 
when the Universe has already flattened out sufficiently 
for mass to be approximately conserved. To define the 
situation more completely, we make two additional mi- 
crophysical assumptions. Firstly, let mass be conserved 
also for individual quantum particles, so that their en- 
ergies scale like E oc c 2 . Secondly, we assume particles' 
wavelengths do not change with c. If homogeneity is pre- 
served, indeed the wavelength is an adiabatic invariant, 
fixed by a set of quantum numbers, eg: A = L/n for a 
particle in a box of size L. 

Under the first of these assumptions a Planckian distri- 
bution with temperature T remains Planckian, but T oc 
c 2 . Under the second assumption, we have A = 2irhc/E, 
and so h/c should remain constant. Therefore the phase 
space structure is changed so that, without particle pro- 
duction, one still has pc 2 cx (fcf,T) 4 /(?ic) 3 , with T oc c 2 . 
A black body therefore remains a black body, with a 
temperature T cx c 2 . If we combine this effect with ex- 
pansion, with the aid of Eqn. (^]) we have 



T + T 







(24) 



We can then integrate this equation through the epoch 
when c is changing to find the temperature Tj of the 
initial state. This fully fixes the initial conditions for 
scalar fluctuations, by means of (p3|). 

In the case of a sudden phase transition we have T + = 



„2 + 



/C 2 



and so 



M 



Mc 2 ~ 



Mc 2 + 



A-(fc) 2 



4fc h T^ 



p^c 



2 + 



(25) 



(26) 



but since A oc c we have 



A+(fc) 



Uk h T+ c+ 



(27) 



Even if T+ = TjS = 10 19 Gew these fluctuations would 
still be negligible nowadays. Therefore although the Uni- 
verse ends up in a thermal state after the phase transi- 
tion, its thermal fluctuations, associated with the canon- 
ical ensemble, are strongly suppressed. 

For a more general c(t) function the procedure is as 
follows. Integrate Eqn. (24) backwards up to a time ij 
when c = 0, to find T(tj). Give A(tj) a thermal spectrum 
of fluctuations, according to j2^), with T(tj). With this 
initial condition integrate Eqns.(|l9|) and (EG) (or even 
better the second order equation fl64| ) given in Appendix 
II), to find A nowadays. 

It is conceivable that a careful design of c(t) would 
leave fluctuations, once c — again, with the right am- 
plitude and spectrum to explain structure formation. In 
particular c(t) may be designed so as to convert a white 
noise spectrum into a scale-invariant spectrum. However 
we feel that until a mechanism for inducing c(t) is found 
such efforts are bound to look ludicrously contrived. 

We feel that the power of VSL scenarios is precisely 
in leaving the Universe very homogenous, after c has 
stopped changing. This would then set the stage for 
causal mechanisms of structure formation to do their job 

sa- 



vin. THE ISOTROPY OF THE UNIVERSE 

There is a sense in which there is an isotropy problem 
in the SBB model, similar to the homogeneity problem. 
We follow closely the remark made in ]13[] , pp.26. 

In Appendix III we write down the vector Einstein's 
equations in the vector gauge, and from them we derive 
the vorticity "conservation" equation when c/c^0. If v 
is the vorticity (defined in Appendix) and II T the vector 
stress, we have: 



a' d 
v' + (l- 3w)—v - 2-v 
a c 



2 1 



-It 1 



w 



(28) 



In the absence of driving stress, v remains constant dur- 
ing the radiation dominated epoch, and decays like 1/a 
in the matter epoch. In |D| it is further argued that the 
relevant dimensionless quantity is 



(k/a)v 
(a'/ca) 



1 



,(l-9™)/2 



(29) 



Hence for w > 1/9 vorticity grows, leading to a further 
fine tuning problem. 

This is most notably a problem if we accept the Planck 
equipartition proposal, introduced in jlT]]. At Planck 
epoch there would then be a significant vorticity. De- 
pending on how one looks at it, this vorticity would then 
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get frozen in or grow, leading to a very anisotropic Uni- 
verse nowadays. 

Whether or not this is a problem is clearly debatable. 
In any case either inflation or VSL models could solve 
this prospective problem. For w < —1/3 we have that v 
decays faster than 1/a 2 . Whatever dimensionless quan- 
tity one chooses to look at, vorticity is therefore safely 
inflated away. If c/c ^ we have that v (x c 2 . Again 
any primordial vorticity is safely suppressed after a phase 
transition in c satisfying any of the conditions (0) or (18). 



IX. THE ENTROPY PROBLEM AND SETTING 
THE INITIAL CONDITIONS 

Let us first consider the SBB model. Let Sh be the en- 
tropy inside the horizon, and ah = Sh/ks be its dimen- 
sionless counterpart. a% is of order 10 96 nowadays. If we 
assume that the only scales in the cosmological model are 
the ones provided by the fundamental constants, then at 
tp the temperature is Tp. At Planck time, ah (being 
dimensionless) is naturally of order 1. In the SBB model 
the horizon distance is dh — It in the radiation domi- 
nated epoch, and ignoring mass thresholds t oc 1/T 2 . If 
evolution is adiabatic one then has (in a flat Universe) 



'T P 

07, (i) ~ O h (tp)[ — 



(30) 



Since anitp) ~ 1, one has <Jh(to) ~ 10 96 . Thus the large 
entropy inside the horizon nowadays is a reflection of the 
lack of scales beyond the ones provided by the funda- 
mental constants, the fact that the horizon size is much 
larger nowadays than at Planck time, and the flatness of 
the Universe. One may rephrase the horizon and flatness 
problems in terms of entropy Q . However if one is willing 
to accept the horizon structure and flatness of the Uni- 
verse simply as features of the initial conditions (rather 
than problems), there is no additional entropy problem. 

There is a problem that arises if one tries to solve the 
horizon problem, keeping the adiabatic assumption, by 
means of superluminal expansion. This blows what at 
Planck time is a region much smaller than the Planck 
size into a comoving region containing the whole observ- 
able Universe nowadays. This solves the horizon prob- 
lem. However if evolution is adiabatic such a process im- 
plies that Chita) <C 1- Stated in another way, since the 
number of particles inside the horizon rih is of the same 
order as ah, this implies an empty Universe nowadays. 

More mathematically, if dh is the horizon proper dis- 
tance, one has 



ah 



3 

dh 



horizon grows faster than t if p < —pc 2 /3: it grows ex- 
ponentially if p = —pc 2 , and like t n (with n > 1) for 
— pc 2 < p < —pc 2 /3. This provides the inflationary so- 
lution to the horizon problem. However in the latter 
case Eqn. d3l| ) implies that ah decreases exponentially, 
leading to OhJto) <C 1. The way inflation bypasses this 
problem is by dropping the adiabatic assumption. Indeed 
during inflation the Universe supercools, and a period of 
reheating follows the end of inflation^]. 

In a VSL scenarios the detailed solution to the entropy 
problem depends on when and what type of "natural con- 
ditions" are given to the pre transition Universe. We first 
derive equations for the entropy under varying c. From 
s = (4/3)pc 2 /T, p (x T 4 /(hc) 3 , and from Eqns § and |l| 
we obtain that the entropy of radiation satisfies 



s 3 p 3 c e(l + £a) 

s 4 p a 2 c 1 + e 



3c 
2c £A 



(32) 



If the Universe is EDSU, there are no violations of mass 
conservation, and entropy is conserved. However if the 
Universe is open or has a positive cosmological constant, 
then we have seen that there is creation of mass. Accord- 
ingly there must be creation of particles, and entropy is 
produced. If the Universe is closed, particles are taken 
away, and the entropy decreases. 

The most suspicious case is therefore if the Universe 
was Einstein de-Sitter before the phase transition. Let 
us assume therefore that at t = tp (the Planck time with 
the constants before the transition) the entropy inside 
the horizon (which has proper size c~tp) was of order 1. 
Then the entropy inside the Hubble volume at t = tp, 
before and after the transition, is 



ah{t + p)=ah{t-p)\ C+t ^ (a{tpy " 



CTtr 



o(t+) 



(33) 



where we have used tp/tp = (c~/c+) 2 . One takes a frac- 
tion (c + /c~) 3 of the horizon volume before the transition 
to make the Hubble volume after the transition. However 
the entropy inside the horizon has increased since tp by 
the same factor. Therefore entropy conservation in this 
case does not conflict with ah(t~^) « 1 after the transi- 
tion. One way of understanding this is that by imposing 
flatness from the outset (before the transition) one has 
already "solved" the entropy problem. Notice that the 
above argument works for any value of t c /tp. 

Now consider the case where "natural" initial condi- 
tions were also imposed at tp, with A = 0. One should 
have e(tp) of order 1. We have already discussed how 
the flatness problem is solved in this case, when large 
empty curvature dominated universes are filled with a 



(31) 



where we have used dh = a f dt'/a. With any standard 
matter (p > —pc 2 /Z) the horizon grows like t. Accord- 
ingly ah grows like a power of t. On the other hand the 



2 This issue has been carefully analyzed in the context of 
inflationary models and models with time varying G in 
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(nearly perfectly) critical energy density during the tran- 
sition. Open Universes become very empty, but they are 
still pushed to EDSU at the transition. One may in- 
tegrate © to find that s+/s~ = (1 + er)- 3 / 4 . One 
may also use Eqn. [To] to find that e has evolved since 
tp to e-(t+) + 1 w (a(t p )/a(t+)) 2 w (t p /t+) 2 , where 
we have used a oc i for the Milne Universe. Hence we 
have that during the transition entropy is produced like 
s+/s- = (tp/t P ) 3/2 = (c" /c+) 3 . Given that a oc f for 
such Universes, the entropy before the transition in the 
proper volume of size c + tp is 

that is there is practically no entropy in relevant volume 
before the transition. However we have that after the 
transition 

a h (t+) = S + (c+tp = S-(c+t+)(^j si (35) 

In such scenarios the Universe is rather cold and empty 
before the transition. However the transition itself re- 
heats the Universe. Notice that, like in the first case 
discussed, the above argument works for any value of 

tc/tp- 

If at t = tp one also has e\ ~ 1 then we have a scenario 
in which the cosmological constant dominates, solves the 
flatness problem, and is discharged into normal matter. 
However if p\ sa pp at t rs tp, then whatever the tran- 
sition time, after the transition the Universe will have a 
density in normal matter equal to p m = pp. Hence the 
Hubble time after the transition will be tp, whatever the 
actual age of the Universe. One may integrate (^) to find 
that in this case (setting e — 0) the entropy production 
during the transition is s + /s~ = (1 + e^) 3 ^ 4 . In the pe- 
riod between t = tp and the transition, e\ increases like 
a 4 , and the entropy density is diluted like 1/a 3 . Hence 
after the transition the entropy density is what it was at 
t = tp, that is s + ss 1/Lp 3 . If we now follow the Uni- 
verse until its Hubble time is ip (when its density is pp) 
we must wait until the expansion factor has increased by 
a factor of (p~p/ Pp) 1 ^ 4 ■ Given that s oc 1/a 3 the entropy 
density is diluted by a factor of (pp/Pp) 3 ^ 4 - Therefore 
the entropy density when the Hubble time is t = tp is 
s sw 1/Lp 3 . Again the dimensionless entropy inside the 
Hubble volume, when this has size Lp~, is of order 1. 

Finally it is worth noting that treating the pre- 
transition universe simply as a Roberston- Walker model 
is no doubt overly simplistic, and we use it simply as 
a device to introducing our ideas. We expect that fur- 
ther development of these ideas could result in a radically 
different view of the pre-transition phase (much as has 
happened with the inflationary scenario). One interest- 
ing observation is that one could avoid having multiple 
Planck times by considering that G oc c 4 . Such assump- 



tion would not conflict with the dynamics of flatness and 
A, as shown before, but now tp = tp. 

X. CONCLUSIONS 

We have shown how a time varying speed of light could 
provide a resolution to the well known cosmological puz- 
zles. These "VSL" models could provide an alternative 
to the standard Inflationary picture, and furthermore re- 
solve the classical cosmological constant puzzle. At a 
technical level, the proposed VSL picture is not nearly 
as well developed as the inflationary one, and one pur- 
pose of this article is to stimulate further work on the 
unresolved technical issues. We are not trying to take 
an "anti-inflation" stand, but we do strongly feel that 
broadening the range of possible models of the very early 
Universe would be very healthy for the field of cosmology, 
and would ultimately allow us to state in more concrete 
terms the extent to which one model is preferred. 

On a more fundamental level we hope to expand the 
phenomenological approach presented in this paper into a 
theory where the concept of (Poincare) symmetry break- 
ing provides the physical basis for VSL. Symmetry break- 
ing is also the central ingredient in causal theories of 
structure formation. We therefore hope to arrive at a 
scenario where symmetry breaking provides a complete 
and consistent complement to the SBB model which can 
resolve the standard puzzles as well as explain the origin 
of cosmic structure. 

Note added in proof: 

After this paper, as well as its sequel f26|| , were com- 
pleted J. Moffat brought to our attention two papers in 
which he proposes a similar idea J27]]. While Moffat's 
work does not go as far as ours in addressing the flat- 
ness, cosmological constant, and entropy problems, he 
does go considerably further than we have in terms of 
specific model building. Moffat's model does not satisfy 
our prescription for solving the cosmological problems, 
but it may do so in a modified form. We are currently 
investigating this possibility. We regret that because we 
were unaware of this work we did not cite it in the first 
publicly distributed version of this paper. 
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APPENDIX I: A SPECIFIC REALIZATION OF 
VSL 

In this Appendix we set up a specific VSL theory. We 
first discuss the simple case of the electrodynamics of the 
point particle in Minkowski space time. We start from 
Bekenstein's theory of variable a, and show how a VSL 
alternative could be set up. We highlight the subtleties 
encountered in the VSL formulation. We then perform 
the same exercise with the Einstein-Hilbert action. We 
briefly consider the dynamics of the field ip — c 4 . Finally 
we cast the key elements of our construction into a body 
of axioms. 



A. Electrodynamics in flat space time 

A changing a theory was proposed by Bekenstein p3| ] 
based on the postulate of Lorentz invariance. The elec- 
trodynamics of a point particle was first analyzed. If 
Lorentz invariance is to be preserved then the particle 
mass m and its charge e must be variable. In order to 
preserve "minimal coupling" (reduction to standard elec- 
tromagnetism when a — const) one chooses the world line 
action 



(36) 



with u M = ± M , g^ v — i]^, e = e(x M ), and m = m{x 11 ). 
Minimal coupling means simply to take the standard ac- 
tion and replace e and m by variables without breaking 
Lorentz invariance. e and m must then be scalar func- 
tions. This action leads to equation: 

{mXfj) = -m^c 2 + ^-u v F^ v (37) 

with the electromagnetic field tensor defined as 

F MV = ^{8 lt {eA u )-d v {eA ll )) (38) 

The electromagnetic action can therefore be defined as: 
-1 



Sem = I d 4 xF^F^ 



(39) 



Also, the particle action ( p6[ ) may be written as a La- 
grangian density: 



Sm 



d^x 



(5( 3 )(x-x(r)) 



7 



(-mc 2 + (e/cju^Af,) (40) 



in which 7 is the Lorentz factor. Maxwell's equations are 
then: 



(41) 



with the current 



3 



x — x(t)) eu^ 



7 



(42) 



This current^] is the current which couples to the gauge 
field, and in the rest frame it equals e. Therefore it can- 
not be conserved, and indeed we have that 



(43) 



Let us now postulate instead that a changing a is to 
be interpreted as c oc Ti oc a -1 / 2 , and that e and m are 
to be seen as constants. Minimal coupling, in the above 
sense, would then prompt us to consider the action (|3^) , 
but with c = c(cc^) everywhere, and e and m constants. 
This action leads to equations: 



1 e 
mx M = -(mc 2 )^ + -u v F^ v 

with the electromagnetic tensor defined as 



F„ 



c(d^A v /c)-d v (AJc)) 



(44) 



(45) 



However the above construction is not complete. In 
spite of the appearance of Eqns. (|36]), (Q), and (|4q) , 
Lorentz invariance is broken. This boils down to the 
fact that, say is no long a 4-vector. Even if c were 
to be regarded as a scalar, would contain c in its 
zero component, but not in its spatial components. The 
usual contractions leading to S could still be taken but 
S would no longer be a scalar. This manifests itself in 
the equations (jl4|) in the fact that in x^ there are terms 
in dc which break Lorentz invariance. 

Since the action is not Lorentz invariant, a minimal 
coupling prescription cannot possibly be true in every 
coordinate system. Minimal coupling is now the state- 
ment that there is a preferred reference frame in which 
the action is to be obtained from the standard action 
simply by replacing c with a field. Let us call this frame 
the "light frame" . In regions in which c changes very lit- 
tle changes in the action upon Lorentz transformations 
are negligible. Hence all boosts performed upon the light 
frame become nearly equivalent and Lorentz invariance 
is recovered. 

The Maxwell equations in a VSL theory become 



-d^(cF^) =47TJ" 



(46) 



in the light frame. Given that Lorentz invariance is bro- 
ken, one can no longer expect the general expression for a 
conserved current to take the form d^j^ — 0. Indeed one 



3 There is an alternative view in which rather than a chang- 
ing e one considers that the vacuum is a dielectric medium 
with variable e. One may then identify a conserved charge, 
but this is not the charge which couples to the gauge field. 
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could try and compute d v of equations ( |46| ) , but now 
and d v do not commute. Also their commutator is not 
Lorentz invariant: for instance [<3o)^i] = (—dic/c 2 )do. 
Still, d^jv = holds in the time frame. It is just that 
this expression transforms into something more compli- 
cated in other frames. The more complicated expression 
would still place constraints on the theory, which could 
still be called "conservation of charge" . 



B. Minimal coupling to gravity 

Let us now examine gravity in such a theory^. As 
in the previous case we will impose a minimal coupling 
principle. Working in analogy with Brans-Dicke theory, 
let us define a field ip — c 4 , and introduce the following 
action 



dx 4 



-/ iP(R + 2A) 
~ g { 16ttG 



C 



M 



(47) 



The dynamical variables are a metric g^ v , any matter 
field variables contained in Cm, and "0 itself. The Rie- 
mann tensor (and the Ricci scalar) is to be computed 
from g M „ at constant ip in the usual way. 

As in the previous section covariance is broken, in spite 
of all appearances, ip does not appear in coordinate 
transformations of the metric, and so the connection r™ 
does not contain terms in Vip in any frame. However the 
connection will contain different terms in ip in different 
frames. Hence the statement that the Riemann tensor 
is to be computed from the metric at constant ip can 
only be true in one preferred frame. Minimal coupling 
requires the definition of a light frame. The action ( |47j ) 
is only Lorentz invariant in appearance. 

Varying the action with respect to the metric leads to: 



6S 
SS M 



8nG 



(48) 
(49) 



leading to a set of Einstein's equations without any extra 
terms 



r a - ^Lt 



(50) 



valid in the light frame. This is the way we chose to 
phrase our postulates in Sectionlll of our paper. In other 



4 Gravitation is normally regarded as the gauge theory of 
the Poincare group |Q. Here we simply abandon this point 
of view. In some future work we will try to define a gauge 
principle for broken symmetries, thereby recovering the stan- 
dard view 



words all we need is minimal coupling at the level of 
Einstein's equations. 

The fact that a favoured set of coordinates is picked 
by our action principle is not surprising as Lorentz in- 
variance is broken. On the other hand notice that the 
dielectric vacuum of Bekenstein theory is an ether the- 
ory. His theory also breaks Lorentz invariance, not at 
the level of the laws of physics, but in the form of the 
contents of space-time. 

In changing a theories a favoured frame is always 
picked up. In a cosmological setting it makes sense to 
identify this frame with the cosmological frame. Free 
falling observers comoving with the cosmological flow de- 
fine a proper time and a set of spatial coordinates, to be 
identified with the light frame. In this frame the Einstein- 
Hilbert action is minimally coupled to a changing ip, and 
the same happens to Friedmann equations. The rest of 
our paper follows. 



C. The dynamics of ip 

The definition of controls the dynamics of ip. This 
is the most speculative aspect of our theory, but it also 
opens the doors to empirical model building. In our 
paper we preferred a scenario in which c changes in 
an abrupt phase transition, but one could also imagine 
c oc a n . The latter scenario would result from a Brans 
Dicke type of Lagrangian 



16nGip 



4' 2 



(51) 



(where us is a dimcnsionless coupling) and is being inves- 
tigated. Addition of a temperature dependent potential 
V(ip) would induce a phase transition, as in the scenario 
developped in our paper. 

However here we only make the following remarks, 
which are independent of any concrete choice of C^p. If 
K = A = one has 



6ip 



(52) 



and so in the radiation dominated epoch (T = 0), once 
K = A = 0, one should not expect driving terms for 
the ip equation. Hence once the cosmological problems 
are solved, in the radiation epoch, c and h should be 
constants. Incidently, once the matter dominated epoch 
(T 0) is reached, ip should perhaps start changing 
again, with interesting observations consequences P2| . 
We are studying the phase space portraits of these cos- 
mologies, when say A ^ 0, and with various C^,. 

During phase transitions the perfect fluid approxima- 
tion must break down. One should then use, say, scalar 
field theory (let's call it <p). Now notice that terms in <p 
will act as a source to ip (as they contain the speed of 
light). Hence whenever there is a phase transition and 
the VEV of a field changes a large amount, one may 
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expect a large change in the speed of light, with most 
choices of C^. A changing ip associated with SSB could 
then solve the quantum version of the cosmological con- 
stant problem, but this might require a rather contorted 
choice of . 



D. Axiomatic formulation of VSL theories 

Postulate 1. A changing a is to be interpreted as 
a changing c and h~ in the ratios c oc Ti oc aT 1 ! 2 . The 
coupling e is constant. 

This postulate merely sets up the theoretical interpre- 
tation of the possible experimental fact that a changes, in 
terms of variable dimensional quantities. This is a matter 
of convention and not experiment, as much as a constant 
Tic is a matter of convention. With the above choice a 
system of units for mass, length, time, and temperature 
is unambiguously defined. 

Postulate 2. There is a preferred frame for the laws 
of physics. This preferred frame is normally suggested by 
the symmetries of the problem, or by a criterium such as 
c = c(t). 

If c is variable, Lorentz invariance must be broken. 
Even if one writes Lorentz invariant looking expressions 
these do not transform covariantly. [In general this boils 
down to the explicit presence of c in the operator d^. 
Once one admits that Lorentz invariance must be explic- 
itly broken then a preferred frame must exist to formu- 
late the laws of physics. These laws are not invariant 
under frame transformation, and one may expect that a 
preferred frame exists where these laws simplify. 

Postulate 3. In the preferred frame one may obtain 
the laws of physics simply by replacing c in the standard 
(Lorentz invariant) action, wherever it occurs, by a field 
c = c(x^). 

This is the principle of minimal coupling. Because the 
laws of physics cannot be Lorentz invariant it will not 
hold in every frame. Hence the application of this postu- 
late depends crucially on the previous postulate supply- 
ing us with a favoured frame. This principle may apply in 
Minkowski space time electrodynamics, scalar field the- 
ory, etc, in which case the frame in which c = c(t) is 
probably the best choice. The cosmological frame, en- 
dowed with the cosmic proper time is probably the best 
choice in a cosmological setting. 

Postulate 4 The dynamics of c must be determined by 
an action principle deriving from adding an extra term 
to the Lagrangian which is a function of c only. 

This is work in progress. We do not wish to specify 
this postulate further because for all we know this extra 
term can be anything. We merely specify that no fields 
(including the metric) must be present in this extra term 
because we wish minimal coupling to propagate into the 
Einstein's equations. 



APPENDIX II: SCALAR PERTURBATION 
EQUATIONS FOR VSL MODELS 

In this Appendix we derive the scalar cosmological 
perturbation equations in VSL scenarios. We assume 
K = A = 0, and use a gauge where the perturbed metric 
is written as 

ds 2 = a 2 [-(l + 2AY)dn 2 - 2BYk i dx i dr) + ,<i.r' <l.r> 

(53) 

for a Fourier component with wave vector k % . Here Y is a 
scalar harmonic. We shall use conformal time rj to study 
fluctuations, and denote ' = d/dr]. The stress energy 
tensor is also writen as 



JT ° = - P Y5 



5TX 



c A J c 



5T] = pU L Y5) + (k l kj - l/3<^r)Y>n T 
The Einstein's constraint equations then read |l| 



ka' 



A 




AirGa 2 

—P 5 

c l 



(54) 



(55) 



B iirGa 2 / p \ v 

C+t - (56) 



and the dynamical equations are 



A+ ±. B' + 2-B) =- 




8irGa 2 p U T 
c 2 c 2 k 2 



(57) 



A 4-KGa 2 p 



(n L -2n T ) (58) 



We assume that these equations do not receive correc- 
tions in c/c. This statement is gauge-dependent, much 
like its counterpart for the unperturbed Eintein's equa- 
tions. We can only hope that the physical result does 
not change qualitatively from gauge to gauge. Comply- 
ing with tradition we now define the comoving density 
contrast |l| 



a! 1 (v 



S + 3(l + w) — -( 

ca k \ 



B 



We also introduce the entropy production rate 

c 2 

r = n L - -^s 

wc z 

where the speed of sound c s is given by 



2 P 2li 

c c = — = wc |1 



2 1 d a 
31 + w c a' 



(59) 



(60) 



(61) 



Note that the thermodynamical speed of sound is given 
by c 2 = (dp/dp)\s- Since in SBB models evolution is 
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isentropic c 2 = (dp/dp)\s = p/p = p'/p'- When c ^ 
evolution need not be isentropic. However we keep the 
definition c 2 = p'/p' since this is the definition used in 
perturbative calculations. One must however remember 
that the speed of sound given in ( |6l| ) is not the usual 
thermodynamical quantity. With this definition one has 
for adiabatic perturbations Sp/5p = p'/p', that is the 
ratio between pressure and density fluctuations mimics 
the ratio of their background rates of change. 

Combining all four Einstein's equations we can then 
obtain the (non-)conservation equations MM 



A' 



a' 

'— 4 




A = 


a 






a' 




).- 


a 






kc 2 


w r 

w 


— kc 



(1 + w)kv — 2 — wTIt 
a 

' elk 



(62) 



1 + w 
2/3 



1 



3 a' 




2k~a 




3 ( 




k 2 c 2 \ 





\wIL T (63) 



These can then be combined into a second order equation 
for A. If r = 11^ = this equation takes the form 



A" + /A' + (.g + / l + c 2fc 2 )A = 



with 



/ 
9 
h 



(l-3w)-3- 
c 

a'V/9 , 



a 



2 W 



3w-- 



9w 
~2~ 



c a 
c a 



(64) 

(65) 
(66) 
(67) 



APPENDIX III: VECTOR PERTURBATION 
EQUATIONS FOR VSL MODELS 

In a similar fashion we can study vector modes in a 
gauge where the metric may be written as 



ds 2 = a 2 [-dr] 2 + 2BY i dx i dr) + Sijdx'dx*] 



(68) 



where Yi is a vector harmonic. The stress energy tensor 
is written as 



5T l3 = P U T Y {1 d) 
Einstein's equations then read p3[ 

16ttG 



k 2 B = 



B' 



2 a -B 
a 



8 -^la 2 Tl T 



(69) 
(70) 

(71) 
(72) 



We assume that these do not receive c/c corrections. The 
conservation equation is then: 
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v + (1 — 3w) — v 
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2—v = 

c 



kc w 
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